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ABSTRACT 

Context. Galactic foreground emission fluctuations are a limiting factor for precise cosmic microwave background (CMB) anisotropy 
measurements. 

Aims. We intend to improve current estimates of the influence of Galactic synchrotron emission on the analysis of CMB anisotropics. 
Methods. We perform an angular power spectrum analysis (APS) of all-sky total intensity maps at 408 MHz and 1420 MHz, which 
are dominated by synchrotron emission out of the Galactic plane. We subtract the brighter sources from the maps, which turns out to 
be essential for the results obtained. We study the APS as a function of Galactic latitude by considering various cuts and as a function 
of sky position by dividing the sky into patches of ~ 15° x 15° in size. 

Results. The APS of the Galactic radio diffuse synchrotron emission is best fitted by a power law, Q ~ k£", with a e [-3.0, -2.6], 
where the lower values of a typically correspond to the higher latitudes. Nevertheless, the analysis of the patches reveals that strong 
local variations exist. A good correlation is found between the APS normalized amplitude, ^loo = kx 100^, at 408 MHz and 1420 MHz. 
The mean APS for £ G [20, 40] is used to determine the mean spectral index between 408 MHz and 1420 MHz, /5(o.408-i.4)ghz ^ 
[-3.2, -2.9] (Q(y) oc y~^^), which is then adopted to extrapolate the synchrotron APS results to the microwave range. 
Conclusions. We use the 408 MHz and 1420 MHz APS results to predict the Galactic synchrotron emission fluctuations at frequencies 
above 20 GHz. A simple extrapolation to 23 GHz of the synchrotron emission APS found at these radio frequencies does not explain 
all the power in the WMAP synchrotron component even at middle/high Galactic latitudes. This suggests a significant microwave 
contribution (of about 50% of the signal) by other components such as free-free or spinning dust emission. The comparison between 
the extrapolated synchrotron APS and the CMB APS shows that a mask excluding the regions with \bgai\ < 5° would reduce the 
foreground fluctuations to about half of the cosmological ones at 70 GHz even at the lowest multipoles. The main implications of our 
analysis for the cosmological exploitation of microwave temperature anisotropy maps are discussed. 

Key words. Galaxy: general - Radiation mechanisms: non-thermal - Methods: statistical - Cosmic microwave background 

1, Introduction forthcoming Planck mission (|Taubed |2QQ4|) . IVloreover, it 

, , . , , n . would enable a precise reconstruction of the individual fore- 

The possibility to study the primordial phases of our Universe ^^^^^ components and therefore a complete astrophysical ex- 
and Its properties and evolution through CMB anisotropies re- imitation of the satellite data. 

lies on our capability to precisely extract the cosmological sig- ^ i ^. , ^ • • • • r 

. r 1 frj. 1 7 1 Ihnnnx a/t ? • Galactic synchrotron emission is the maior source or con- 

nal from observations (Tegma rk et al.n 2000). JVlaps or the mi- ^ . . ^-V • i i rr. ^tt r • . j- . j 

1 . 1 J > 1 \i . 1 .1 tammation at frequencies below 50 GHz for intermediate and 

crowave sky include many Galactic and extragalactic astrophys- , , ■, i . 

. 1 ^ ^^: ^- A ^ r r^K/TT^ • ^ large angular scales, as recently confirmed by the impres- 

ical contributions. A correct recovery of the CMB anisotropy . ^ ...nv/TAo rB TTTi rrr — u — ^ 1 1 hnn-y K 

^ 1 , . ^ 1 r r A ' ^ sive WMAP results (Bennett et al. 2003, Hmsha w et al.ll2QQ7D . 

field requires an accurate removal or those foreground signals ^ , ^ i-. /t^i i ot-i. ^r^r^r^. — — 

r .11. 1 .1 ij rSr j Synchrotron radiation Rvbicki & Lishtman 1979) ariscs from 

from the observed maps. Current knowledge of the foreground . , ^ ' • • ..cur 

^ ^ ^ . .1 u 11 r cosmic ray electrons gyrating m the magnetic field of our 

components permits one to retrieve the bulk of the cosmo- ^ , ^/ .1 f r .1 

, ' , • r ^- J J • r^A/TT^ • ^ J • Galaxy. The energy spectrum and the density of the cosmic ray 

logical information encoded m the CMB anisotropy, and, m , ^ n u . .1 .1 

^,•1 1 , / A T^ox TVT .1 1 electrons as well as the magnetic field strength vary across the 

particular, its angular power spectrum (APS). Nevertheless, a ^ , ^, ^ ^, 11 1 . • • -n j 

J J.J- r.u r J • • 1 . ..1 • Galaxy, therefore the observed synchrotron emission will de- 

deeper understanding of the foregrounds is crucial to settle im- 1 ' 1 ^ , .1 • r .1 1 j • 1 

5 ^ I I- u • u r TT^A/rAr»m pend on the frequency and on the region of the sky. Radio ob- 

portant cosmological issues, which arose from the WMAP M ^ . ^ r V^tt j .7 1 . • . r .1 

^ 1. 11.1 onn/c XT 11 o 1 u J ^nn/T scrvations at y < 5 GHz provide the clearest picture of the 

results (Naselskv et al. 200 6; Naselskv & Verkhodanov 2006; ^ , ,7 ,1 . . r . 

IT^TT — iiionA-zK J ij . 11 u J J J u Galactic synchrotron morphology, since at these frequencies the 
IChiang et al.l I2QQ7D and could potentially be addressed by the , ^. 11, . . n .1 

^ r- J J diffuse non-thermal radiation clearly dominates over all other 

Send offprint requests to: L. La Porta 

^ http://lambda.gsfc.nasa.gov/ ^ http://www.rssd.esa.int/Planck 
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emission components outs ide the Galactic pla ne. In the past, 
the 408 MHz all-sky map fHaslam et al.l ll982h has extensively 
been used as a template for the Galactic synchrotron emission in 
foreground separation attempts (e.g. Bouchet & Gispert 1999; 
Bennett et al. 2003). In addition, that map as well as other 
less suited surveys have been exploited to find an appropriate 
parametrization of the synchrotron emission APS to be /) extrap- 
olated to the microwave range for estimating the contamination 
of CMB anisotropics at diff'er ent angular scales dGiardino et al.l 
I2OOII: iBaccigalupi et al.|[200ll) and //) used to se t priors in fore- 
ground separation applications ( Tegmark & Ef stathioiJ Il996l: 
iBouchet et anil999l: iBouchet & Gispert 1999). The outcome of 
these analyses is that the APS of the synchrotron emission com- 
puted over large portions of the sky can be modelled by a 
power law, i.e. Q ~ ^ ~ 180°/^), with a spectral index 

a ~ [-3, -2.5] for ^ < 200, corresponding to angular scales 
^ > 1°. Clearly, such a general result is just a first step as it does 
not describe the complexity of the synchrotron emission APS, 
whose parameters are expected to change with frequency and 
with sky direction. 

We have carried out a detailed analysis of all-sky radio maps 
and improved on previous attempts by providing an accurate 
characterization of the synchrotron emission APS. The results 
obtained for the new 1.4 GHz polarization all-sky survey (Reich 
et al., in prep.) will be reported in a forthcoming companion 
paper. The analysis presented in this paper focuses on the 
synchrotron emission APS in total intensity. A substantial im- 
provement was possible by using a new all-sky map at 1 .42 GHz 
(Reich et al., in prep.), which has a higher angular resolution 
and better sensitivity, in addition to the all-sky map at 408 MHz. 
These maps are currently the best suited data for studying the 
Galactic synchrotron emission at largest angular scales. A more 
detailed description of some technical aspects related to this 
work can be found in iLa Portal ( 2007h. W e extensively used 
the HEALPixEl software package (iGorski et al. 2005). HEALPix 
(Hierarchical Equal Area isoLatitude Pixelization) is a curvilin- 
ear partition of the sphere optimized for fast spherical harmonics 
transforms and angular power spectrum estimation. The latter 
task is performed by the facility Anafast. We also made use 
of the data reduction package based on the NOD2-software 
(lHaslamlll974l) . 



The paper is organized as follows. Sect. 2 describes the char- 
acteristics of the 408 MHz and 1.42 GHz total intensity surveys, 
their projection onto HEALPix maps and noise considerations. 
In Sect. 3 the Galactic radio emission APS over large areas 
is examined, which reveals the necessity of a discrete source 
subtraction for a correct evaluation of the diff'use synchrotron 
APS. The two all-sky maps are decomposed into a map of the 
diff'use component and a map of discrete sources. Their angular 
power spectra are derived and discussed. The results obtained by 
fitting the angular power spectra of the diff'use component maps 
are presented. In Sect. 4 the radio survey angular power spectra 
are extrapolated to the microwave range for a comparison 
with the WMAP 3-yr results. Sect. 5 is dedicated to a local 
analysis of the radio map APS. We summarize our results and 
conclusions in Sect. 6. 



2. The data 

The present analysis focuses on the APS of the Galactic 
synchrotron emission at radio frequencies. Ho wever, the 23 
GHz synchrotron component obtained by Hinshaw " et al.l (l2007h 
using the WMAP 3-yr data has also been considered to some 
extent and will be further discussed in Sect. IH 



2.1. The 408 MHz and 1420 MHz surveys 



The 408 MHz map ('Hasla m et al.lll98 2^ was produced by merg- 
ing diff'erent datasets obtained with large parabolic reflector tele- 
scopes (Jodrell Bank 76 m, Efl'elsberg 100 m and Parkes 64 m 
telescopes - see Fig. 1 of Haslam et al. 1982), using a similar 
observing strategy and the same calibration procedure. The final 
map is characterized by an angular resolution of Ohpbw ~ 0!85 
and a 20' pixel rms-noise of about 670 mK. The version used 
in the present analysis was corrected for a zero level problem 
concerning the portio n of the sky observed from Jodrell Bank 



concerning the portio n 
(iReich & Reichl(l988h . 



The total intesity map at 1420 MHz has been obtained by 
combining a northern and a southern sky survey (Reich et al., 
in preparation). Both surveys are on an absolute temperature 
scale a nd zero level by usi ng low resolution sky horn measure- 
ments (iTestori et al.ll200ll) . This includes a correction for far- 
side lobe contamination for single-dish telescopes. Receiving 
systems were used, which provide total intensities unaff'ected 
by linear polarization. The northern sky survey was observed 
with the Stoc kert 25-m telescope and extends in d eclination from 
-19° to 90° (lReichlll982l: iReich & Reichlll986l) . The southern 
sky survey was carried out with the Villa Elisa 30-m telescope 
i n Argentina and covers the latitude range 6 e [-90° to -10°] 
(IReich et al.l200lh . Both have an angular resolution Ohpbw ~ 36' 
and overlap for latitudes in the range [-19°,- 10°]. The result- 
ing map has a 15' pixel rms-noise of ~ 17mK. The original 
maps are provided as equidistant cylindrical (ECP) maps. For 
the presen t analysis these rn aps have been projected into the 
HEALPix (iGorski et al.l l2005l) pixelization scheme adopted by 
the WMAP and Planck Consortia. For this purpose a simple 
regridding al gorithm has been imp lemented, which is described 
in detail by iLa Porta et al.l (120051) . The reliability of the pro- 
jection provided by this simple approach has been verified by 
successfully performing forward and backward transformations 
between the two tessellation schemes. The produced HEALPix 
maps have a pixel size of ~ 1' (the number of pixels for an all- 
sky map is Npix = 12 here we used riside = 512). 



2.2. Noise estimate 

The authors of the radio maps have estimated the rms-noise di- 
rectly on the ECP maps by calculating the standard deviation 
of low emission regions. Going from a Cartesian representation 
of the sphere to the HEALPix scheme, the rms-noise per pixel 
should theoretically decrease toward the polar caps, according 
to the formula: 



1 



X Cr pixel,ECP 



http://healpix.jpl.nasa.gov/ 



where N is the number of the ECP pixels corresponding to each 
HEALPix pixel at a fixed latitude. Such a relation holds under the 
hypothesis that the noise is Gaussian and uncorrected among the 
pixels. Both these assumptions are not necessarily satisfied in the 
examined surveys, in primis because the pixel size is about half 
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Fig. 1. Angular power spectra of the Galactic plane cut-offs 
{\hgai\ > hcut) at radio frequencies. Each line refers to a certain 
bcut (from the top: black ^ 40°, fuxia ^ 50°, green ^ 60°). 



the angular resolution. Furthermore, the rms-noise quoted for the 
408 MHz and 1420 MHz surveys quantifies the temperature fluc- 
tuations per pixel, which is due not only to the receiver white 
noise, but also to the contribution of unresolved sources and of 
residual systematics efl'ects (as "scanning strategy efl'ects"), that 
could not be fully eliminated in the data reduction procedure. 
The overall rms-noise in the HEALPix maps is therefore proba- 
bly higher than the guess determined in this way, thus making 
this formula suitable for deriving a lower limit. With this for- 
mula we constructed a map of rms-noise at both frequencies. The 
rms-noise per pixel decreases for increasing latitude and varies 
between ~ 10 mK and ~ 18.7 mK at 1420 MHz and ~ 0.5 K 
and ~ 0.7 K at 408 MHz. Due to varying integration times, the 
rms-noise is not constant over the ECP map, but it is expected to 
diminish toward the celestial poles. Also taking this efl'ect into 
account, while applying the above formula, the derived estimate 
decreases at most by a factor ~ 2. This way we obtain an opti- 
mistic lower limit for the rms-noise in the HEALPix maps. The 
corresponding APS is computed as: 



' An (TIN. 



pix 



where cr is the mean value of the rms-noise in the considered 
area and Npix is the total number of pixels in the HEALPix map. 
A generous upper limit for the noise APS is provided instead by 
the high multipole tail of the APS. This way we have bracketed 
the noise APS for each considered sky region. 

2.3. Statistically significant multipole range 

The HEALPix maps produced at 408 MHz and 1420 MHz con- 
tain precise statistical information only for angular scales that 
are: 

- larger than the beam, i.e. for 6 > Ohpbw, therefore the maxi- 
mum multipole relevant in the APS analysis of these surveys 

is ^max ~ 180°/^//pfi^[°]; 

- smaller than the maximum angular extent of the considered 
area, Ocov, so that the coverage sets the minimum multipole. 
A safe choice is ~ 5 x 180°/^cov[°]- 

3. The APS of large areas: analysis for various cuts 

To investigate the dependence of the synchrotron emission APS 
parameters on latitude, the Galactic plane has been "cut ofl'" 
from the maps at difl'erent latitudes by setting to zero pixels with 




10 1 00 1 1 00 1 1 00 1 1 00 

Multipole (/) Multipole (/) Multipole (/) Multipole (/) 

Fig. 2. Comparison between the cut-off angular power spectra of 
the original (fuxia lines at the top) and destriped (blue) version 
of the 408 MHz map. The cut-ofl' angular power spectra of the 
difl'erence map are also shown (green lines at the bottom). 



\bgai\ < bcut, where bcut = 5°, 10°, 20°, 30°, 40°, 50°, 60°. At the 
same time, this approach preserves the largest possible cover- 
age, important to keep the widest range of statistically significant 
multipoles. We also considered "asymmetric cuts", i.e. regions 
with bgai > bcut (northern cuts) and bgai < -bcut (southern cuts), 
thus taking into account the diff'erence between the two Galactic 
hemispheres. In fact, the northern hemisphere contains a large 
and bright feature of the radio sky, i.e. the North Polar Spur 
(NPS), which is widely believed to be an old supernova remnant 
in the Solar System vicinity (gaiter 1983; Egger & Aschenbaclj 

[laH). 

We computed the corresponding APS by using the facility 
Anaf ast of the HEALPix package and renormalized it to account 
for the incomplete sky coverage. The angular power spectra de- 
rived for \bgai\ > bcut with bcut ^ 40° are shown in Fig. [T] as 
representative examples. All the recovered angular power spec- 
tra flatten towards higher multipoles. Such a behaviour of the 
APS might be due to noise, systematic efl'ects ("stripes"), dis- 
crete sources or might be an intrinsic characteristic of the syn- 
chrotron emission fluctuation field. Instrumental white noise can 
be discarded because its APS should be constant, whereas after 
the flattening the angular power spectra decrease with £ as in the 
presence of beam smoothing. "Stripes" are systematic baseline 
distortions in the telescope's scanning direction and are mainly 
due to the limited stability of the receiving system and to the 
influence of weather conditions. They could also be cancelled 
from the list of possible causes from the comparison with a de- 
striped version of the 408 MHz map (Platania et al. 2003). The 
cut-ofl' APS of the two versions of the 408 MHz map present 
only marginal difl'erences at intermediate multipoles (see Fig. [2]). 



3.1. Discrete source subtraction 

Beside difl'use emission, a large number of discrete sources 
(DSs) are visible in the radio maps. A DS subtraction has been 
done by performing a 2-dimensional Gaussian fitting that also 
provides an estimate of the diff'use background, which is approx- 
imated by a tilted plane. Such an estimate has been used to fill 
the pixels originally corresponding to the subtracted DSs. Where 
the background emission shows strong gradients the source fit- 
ting is more difficult. Consequently, the flux limit above which 
all discrete sources most likely have been subtracted is difl'er- 
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TI @ 408 MHz (n^^dg = 512), DSs subtracted map 




10500 ^^^^H ^^^H 80000 (mk) 



TI @ 1420 MHz (113,^^ = 512), DSs subtracted map 




2500 ^^^^^^^^K ^^^^^m 5500 (mk) 



Fig. 3. Mollweide projection of the HEALPix maps produced at 
408 MHz (top) and 1420 MHz (bottom) by subtracting the dis- 
crete sources. The maps are in Galactic coordinates. The center 
of the maps is / = 0, /? = 0. The Galactic longitude increases 
toward the left up to 180°. 



Galactic plane cut — off at b = 5° Galactic plane cut — off at b = 20° 




Fig. 4. HEALPix maps of the Galactic plane cut-offs extracted 
from the 1.4 GHz all-sky map after DS subtraction. A moll- 
weide projection is displayed for the cuts with \bcut\ ^ 20° and 
a gnomonic view (centered on a Galactic pole) for those with 
\bcut\ > 40°. 



ent close to the plane and far out of it. Namely, for \b\ > 45° all 
the sources with peak flux above ~ 0.8 Jy (respec. ~ 6.4 Jy) have 
been subtracted from the 1420 MHz (respec. 408 MHz), whereas 
for I Z? I < 45° such a source detection threshold is ~ 4.6 Jy (re- 
spec. ~ 63.8 Jy). All discrete sources that could be reasonably 
fitted by a Gaussian profile have been eliminated and two new 
maps have been generated at 408 MHz and at 1420 MHz (Fig.O. 

The subtracted DSs (see Burigana et al. 2006 for a map) are 
mostly point sources, except for some rather extended objects, 
as for instance the radiogalaxy Centaurus A, appearing in the 
original radio maps right of the Galactic center at bgai ~ 20°. 
Such extended objects are among the brighter subtracted DSs 
and concentrate along or in the proximity of the Galactic plane. 
They are mainly Galactic sources, i.e. Hll-regions or supernova 
remnants. On the contrary, the DSs subtracted at medium and 
high latitudes are nearly all extragalactic sources. 

The maps of some Galactic plane cut-off's at 1420 MHz (re- 
spec. 408 MHz) after DS subtraction are shown in Fig. [4] (respec. 
Fig. [5]). The maps are displayed adopting diff'erent scales, in or- 
der to emphasize the relative importance of the various com- 
ponents. Note that "scanning strategy eff'ects" are clearly visi- 
ble in the southern sky at high latitude in the map at 408 MHz. 
As discussed above, the angular power spectra of the destriped 
and original versions of the Haslam map do not exhibit signifi- 
cant diff'erences, thus implying that "stripes" are not an issue for 
the APS analysis at 1420 MHz. Figure [6] (respec. Fig.jT]) shows 
the APS of the Galactic plane cut-off's at 1420 MHz (respec. 
408 MHz) for the original map, for the DS-subtracted map and 



for the map of subtracted DSs. The APS of the DS maps almost 
perfectly matches the ffat part of the original map APS at large 
L This result identifies DSs as the reason for the ffattening of the 
original APS and also confirms that the major contribution from 
source contamination has been eliminated in the DS-subtracted 
maps. At high latitude the APS of the Galactic ffuctuation field 
is dominated by the DS contribution for £ > 100, which is due 
to the enhanced relative contribution of the DSs respect to the 
weak diff'use background emission. 

Figure[8] shows the APS of the various cuts for the DS map at 
1420 MHz. Note that for all the southern cuts and for the north- 
ern cuts with bcut ^ 30° the DS angular power spectra are rather 
ffat up to ^ ~ 100 and then decrease as for a beam smoothing. 
On the contrary, for the northern cut at 5° - 20° the DS angu- 
lar power spectra present a power law behaviour at lower mul- 
tipoles, that implies the existence of signiffcant ffuctuations also 
at the larger angular scales, as expected in presence of relatively 
extended discrete structures. The DS APS of the cut at 20° is su- 
perimposed on that of the cut at 30° for £ > 100, whereas it ex- 
hibits a power law behaviour at lower multipoles. The diff'erence 
between the angular power spectra of the northern cuts at 20° 
and 30° is then due to the DSs located in the portion of the sky 
characterized by 20° < bgai < 30°, which includes Centaurus A, 
thus explaining the power law behaviour of the APS at the lower 
multipoles. The same situation was found at 408 MHz. 
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Galactic plane cut — off at b = 5° Galactic plane cut — off at h = 20° 




Fig. 5. As in Fig.H but at 408 MHz. 
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Fig. 6. APS of some Galactic plane cut-off for the 1420 MHz 
maps: (from the top in each panel) original fuxia), after dis- 
crete source subtraction blue) and DSs only green). 



3.2. The APS after source subtraction 

The angular power spectra of the maps after source subtrac- 
tion approximately follow a power law, as expected for the dif- 
fuse Galactic synchrotron emission. For \bcut\ < 40° the angu- 
lar power spectra of all symmetric (\b\ > hcut) and asymmet- 
ric {b < -bcuu b > bcut) Galactic cuts are very similar to each 
other and appear progressively shifted downward (see top panels 
of Fig. O. This result reflects the fact that the Galactic difl'use 
emission becomes weaker for increasing latitude. The angular 
power spectra of the symmetric cuts with \bcut\ ^ 40° are super- 
imposed. The same result holds for the angular power spectra of 
the northern cuts, whereas in the southern hemisphere the APS 



amplitude decreases for increasing \bcutV This discrepancy leads 
to the conclusion that the angular power spectra of the symmet- 
ric cuts with \bcut\ ^ 40° are mainly influenced by the northern 
hemisphere. Indeed, the angular power spectra of the northern 
cuts at bcut ^ 20° have amplitudes larger than those of the south- 
ern cuts at both frequencies. As an example. Fig. [TOl shows the 
comparison between the angular power spectra of the northern 
and the southern cuts at 1420 MHz. The angular power spectra 
of the two Galactic hemispheres can reasonably be expected to 
be similar for the Galactic difl'use synchrotron emission, while 
they turn out to be difl'erent in amplitude and to some extent 
(mostly at smaller scales) in shape. That difl'erence results from 
the combination of two efl'ects. In the southern sky, the angu- 
lar power spectra of the DS-subtracted maps tend to flatten at 
i ~ 150 - 200 due to the presence of unsubtracted sources, 
whose relative contribution to the fluctuation field increases be- 
cause of the low background signal. In the northern hemisphere, 
the Galactic diff'use synchrotron emission is strongly influenced 
by the radiation of the NPS. 

3.3. Fit of the APS after source subtraction 

The radio maps after source subtraction include two astrophysi- 
cal components: the Galactic difl'use emission and the (mainly) 
extragalactic source contribution, which are convolved with the 
telescope beam and contaminated by the instrument noise, that 
can be approximately treated as white noise. We therefore ex- 
press the corresponding APS as 



^map 



/^synch 



)W^-hc"^'^^ (1) 
^(^+1)0-2 ^j^^ window function of the symmet- 



where Wi = 

ric and Gaussian beanfl with cr/^ = Ohpbw[i"<^^]/ V81n2). The 
synchrotron emission APS is empirically modelled as C^^"^^ = 
ki^. We note that such an empirical choice qualitatively can 
be explained by ma^ netohydrod ynamic turbulence arguments 
(IChe p urnov'T998: iCho & Lazarianll2002L 120031) . The contribu- 
tion of the unsubtracted DSs is approximated by a constant term, 
according to the forma Hsm of Poisson fluctuations from extra- 
galactic point sources (iFranceschini et al.lil98 9). The contribu- 
tion of fluctuations due to source clustering is expected to be 
negligible with respect to the Poisso n term at the source detec- 
tion threshold achieved in our maps ( Toff'olattietal.lll998l) . 



^ The pixel window function has been also taken into account, but 
its effect is not important here, because the pixel size is significantly 
smaller than cr^. 




Fig. 7. As in Fig.[6l but at 408 MHz. 
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Fig. 8. Angular power spectra of the northern (bgai > bcut 
left panel) and southern {bgai < -bcut right) cuts for the 
map of discrete sources at 1420 MHz. Color legend (see online 
version): black (dotted) \bcut\ = 5°, black \bcut\ = 10°, 
green \bcut\ = 20°, red \bcut\ = 30°, dark blue \bcut\ = 
40°, fuxia ^ \bcut\ = 50°, light blue ^ \bcut\ = 60°. 
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Fig. 9. Comparison between the APS of the cuts for the DS- 
subtracted map at 1420 MHz. The left (respec. right) panels dis- 
play the angular power spectra of the northern (respec. southern) 
cuts. First (respec. second) row panels: bcut = 5° - 40° (respec. 
i^cut = 40° - 60°) . Color legend (see online version): black dot- 
ted ^ \bcut\ = 5°, black ^ \bcut\ = 10°, green ^ \bcut\ = 20°, red 
^ \bcut\ = 30°, dark blue ^ \bcut\ = 40°, fuxia ^ \bcut\ = 50°, 
dark green \bcut\ = 60°. 



In order to derive the range of variability of the synchrotron 
emission amplitude and slope, two extreme cases have been con- 
sidered (see Appendix C of La Porta 2007 for details). The flat- 
test synchrotron APS compatible with the data is found by ne- 
glecting the source term in Eq. [T] and the steepest one is recov- 
ered by assuming a null noise contribution and maximizing the 
source term. We performed a least-square fit to the APS by ex- 
ploring the parameter space on adaptive grid^. The uncertain- 
ties on the best fit parameters are derived as the diff'erence with 
those obtained in the two extreme cases. Fi gures [TTI and [T2l show 
the angular power spectra and the best fit curves corresponding 
to the best model at the two frequencies, while Tables [T] and [2l 
list the obtained parameters and their uncertainties. For the syn- 



^ For this purpose we implemented a specific algorithm and tested 
its reliability with th e MINUITS package of the CERN libraries 
dJames&RoosIl 19751) . 
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Fig. 10. Comparison between the angular power spectra of the 
northern (bgai > t^cut^ blue - upper lines) and southern (bgai < 

~bcut 

MHz. 



red - lower) cuts for the DS-subtracted map at 1420 



chrotron term, we quote the value of the normalized amplitude 
^100 = ^ X 100^, which corresponds to a physical quantity. In 
fact, ^100 = Q=ioo, thus implying that the normalized ampli- 
tude gives the mean temperature fluctuations at angular scales 
of ~ 2°. Figure fT3] shows the best fit parameters of the syn- 
chrotron APS as a function of the Galactic latitude. The nor- 
malized amplitude, kioo, is maximum when the considered cut 
includes the lower latitudes, where the Galactic radio emission 
peaks. In particular, at 408 MHz (respec. at 1420 MHz) ^loo e 
[488, 6527] mK^ (respec. [0.21, 4.57] mK^) for the northern cuts 
andy^ioo e [138,6734] mK^ (respec. [0.03,4.00] mK^) for the 
southern cuts. The mean error on ^loo is of ~ 18% for the cuts at 
the lower latitude (\bcut\ e [5°, 30°]) and of ~ 30 - 40% for the 
others. The uncertainty is larger for the cuts at higher latitude due 
to the reduced multipole range suitable for the fitting procedure. 
The slope of the synchrotron APS for the northern cuts varies 
in the interval ~ [-3.0,-2.8] at 408 MHz and ~ [-3.0,-2.7] 
at 1420 MHz, while in the southern cuts a ~ [-2.9, -2.6] and 
a ~ [-3.0,-2.7], respectively. The errors on a are on average 
~ (5 - 7)% for \bcut\ < 30° - 40° and typically increase to ~ 18% 
for cuts at higher latitude. 

At both frequencies there is no evidence of a systematic depen- 
dence of the synchrotron emission APS slope on latitude (see 
Fig.m. 
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Fig. 11. Angular power spectra of the northern (top panels) and 
southern cuts (bottom) for the DS-subtracted map at 1420 MHz 
together with the best fit curves obtained. The individual con- 
tributions of synchrotron emission (blue lines) and of sources 
(green lines) are also plotted (smoothed by the beam). 
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Fig. 12. As in ¥ig.\n\ but at 408 MHz. 



The source term c^^^ increases for decreasing latitude, as 
expected given that the source subtraction is less complete be- 
low ~ 45°. F rom extragalactic sou rce counts at 1.4 GHz in 
total intensity (iPrandoni et al.ll200ll) . the expected source con- 
tribution is c^^^ ^ 0.06mK^ for flux densities below ~ 1 Jy 
^^src _ [0.03-0.3] mK^ including the quoted Icr errors and 
also considering the effect of the finite sampling) and c^^^ ^ 
0.30 mK^ (c''' ^ [0.15 - 1.50] mK^ ) for flux den sities be- 
low ~ 5 Jy. At 408 MHz the available source counts (iJamrozyl 
12004 lead to c''' ^ 200 mK^ for flux densities below ~ 6.4 Jy 
(c''' ^ [150 - 360 mK^) and c''' ^ 660 mK^ for flux densities 
below ~ 64 Jy (c''' ^ [260 - 1150] mK^). The values resulting 
from the fits are consistent with the above estimates. 



Table 1. Best fit parameters obtained by modelling the angular 
power spectra of the northern and southern cuts at 1420 MHz 
according to Eq. [T] 
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Fig. 13. Best-fit parameters obtained for the Galactic radio syn- 
chrotron emission APS against Galactic latitude. 



Table 2. As in Table fflbut at 408 MHz. 



4. Extrapolation to the microwave range 

In this section we extrapolate the results obtained from the 
analysis of the 408 MHz and 1420 MHz surveys over large areas 
to the microwave range, in order to make a direct comparison 
with the WMAP 3-yr results. The main objective of the WMAP 
mission was the realization of a CMB anisotropy total intensity 
map and the estimation of the corresponding APS. By-products 
of the mission are maps of the foregrounds contaminating 
the cosmological signal at the five frequencies observed by 
the satellite (v ~ 23,33,41,61,94 GHz), i.e. the microwave 
emission from the Milky Way, characterized by diff'use (dust, 
free-free, synchrotron) and discrete (e.g, HII regions, SNRs) 
components, and from extragalac tic sources. The maps were 
worked out (iHinshaw et al.l l2Q07l) by using templates of the 
various astrophysical components, constructed by exploiting 
ancillary data. Then a pixel-by-pixel (MEM based) fit of all the 
maps (see Bennett et al., 2003 for a description of the method), 
i.e. templates and WMAP frequency maps after subtraction of 
the CMB anisotropy field, was performed posing some priors 



60 



8 



La Porta et aL: APS analysis of total intensity all-sky surveys 



on the spectral behaviour of the foregrounds. 



4.1. Comparison with the WMAP K-band synchrotron 
component 

The 2 3 GHz (K-band) synchrotron map by iHinshav^ et alJ 
((20071) is considered in this section for a comparison with the 
408 MHz and 1420 MHz data. It is evident that such a map 
provides a picture of the global non-thermal emission observed 
by the satellite at 23 GHz, rather than the Galactic synchrotron 
component only. Several extragalactic sourc es are clearly 
recognizable. Furthermore, as pointed out by iHinshaw et al.l 
(l2QQ7h . the diffuse non-thermal emission is concentrated at low 
latitudes and appears remarkably well correlated with the dust 
component 0. This might suggest the presence at 23 GHz of 
anomal ous dust emission, as due for example to spinning dust 
grains (iDraine & LazarianI Il998l) . In fact, most dusty active 
star-forming regions are localized along the Galactic plane. This 
hypothesis seems further supported by the joint analysis of the 
WMAP maps (1-yr release) and the Green Bank Galactic Plane 
Survey by Finkbeiner (2004). De Oliveira-Costa et al. (2004) 
estimated the fluctuations expected at 10 GHz and 15 GHz for 
the foreground component traced by the K-MEM synchrotron map 
by Bennett et al. (2003) (1-yr results), by cross-correlating the 
latter with the Tenerife 10 GHz and 15 GHz CMB maps and all 
the WMAP CMB maps. They found values one order of magni- 
tude below what is expected for the synchrotron emission and 
concluded that the K-MEM synchrotron component by Bennett 
et al. (2003) is dominated by anomalo us dust emission even at 
\bgai\ ^ 20°. iHildebrandt et al.l (l2007h also found evidence for 
anomalous microwave emission at high Galactic latitudes by 
cross-correlating the COSMOSOMAS 11 GHz observations 
with the WMAP K- and Ka-ban d map. The same conclusion 
also was reached by iDavies et"an (l2QQ6h . who cross-correlated 
the WMAP 1-yr map with foreground templates in a dozen 
small patches located at medium and high latitude. However, 
the origin of the spatial correlation found at WMAP frequencies 
between synchrotron and dust emission is still a matter of 
debate. Bennett et al. (2003) claim that the observed correlation 
is the result of a spatially varying synchrotron spectral index, 
which significantly alters the morpholo gy of the synchrotron 
emission with frequency. Hinsha w et al.l (120071) affirm that the 
issue is left open also by the WMAP 3-yr results and that high 
quality and large coverage surveys at v ~ 5 - 15 GHz are needed 
for a decisive test of both the above discussed explanations of 
the synchrotron-dust correlation. 

The separation of the free-free and synchrotron emission in 
low latitude regions is also very uncertain. On one hand, the 
408 MHz map used as a template of the non-thermal emission 
contains a non-negligible contribution (< 10%) of f ree-free at 
lower latitudes (Dickinson et al. 2003; Paladini et al l' 2005). On 
the other hand, the H -a map used as template for the free-free 
emission (Fi nkbeinerl l2003h cannot be properly corrected for 
dust extinction for \bgai\ < 5°, thus potentially leading to an 
underestimation of the expected thermal emission at 23 GHz. 
The situation in the vicinity of the Galactic plane is extremely 
complicated and remains unclear. 

The extrapolation of the angular power spectra derived 
at 408 MHz and 1420 MHz to the microwave range is a 



^ Such a tight synchrotron-dust correlation holds at all WMAP fre- 
quencies. 



delicate issue. The astrophysical components contributing 
to the fluctuation field APS of the radio surveys scale with 
frequency in a diff'erent way. For the Galactic radio emission 
between 408 MHz and 1420 MHz [Reich et al.l (l2004l) compiled 
a map of the spectral index /3 (Tt, oc v^) that reveals a complex 
structure, due to superposition of the spectral behaviour of the 
map components (synchrotron emission, sources, free-free). 
The situation is further complicated by a possible but not well 
known steepening of the diff'use synchrotron emission power 
spectrum above 10 GHz, due to the steepening of the c osmic 
ray electron energy s pectrum (iBanday & Wolfendald Il990l 
1 19911: ISti-ong et al.''2007). Last but not least, the fact that the 
astrophysical components of the map scale with frequency in 
a diff'erent way may imply a change in the overall APS shape, 
since the relative weight of the foreground contribution to the 
APS could vary significantly. 

Given the complexity of the open issues discussed above, 
the following analysis merely aims to verify the consistency 
between the information about the non-thermal radiation APS 
coming from the 408 MHz and 1420 MHz data and from the 23 
GHz WMAP data. We focused on what happens at medium and 
high Galactic latitudes, since the problems in interpreting the 
Galactic emission are more complicated close to the plane. 

We first carried out a source subtraction on the 23 GHz map 
at intermediate and high latitudes (\bgai\ > 40°), similar to that 
performed on the radio surveys. The comparison of the map of 
subtracted sources with the mask of sources produced by the 
WMAP team shows that most (~ 80%) of the objects have been 
identified and subtracted. We have derived and compared the 
APS of the original, source- subtracted and source map at 23 
GHz for some northern and southern cuts. Namely, we consid- 
ered \bcut\ = 40°, 50°, 60° and verified that cf^' """"^ ~ c''' over 
the significant multipole range (i.e. £ > 20). This means that in 
the 23 GHz map the source contribution dominates the high lat- 
itude cut APS at all angular scales. In Figure [141 we show the 
APS of the 23 GHz map after source subtraction for the asym- 
metric cuts with bcut = ±40°. For comparison, we also display 
the extrapolated radio APS (from Tables [T]and[2]), derived as 

Q(23) = QiVradio) • (23/yradiof^ 

where Vradio = 0.408, 1.420 GHz. The frequency spectral index 
j3 is chosen case by case as the value that brings the radio APS to 
overlay the WMAP one at the lower multipoles (the exact values 
are reported in the figure caption). The extrapolated angular 
power spectra are very similar to each other, but steeper than 
those of WMAP. Furthermore, the frequency spectral indices 
needed in the extrapolation suggest the existence of a steeper 
spectral behaviour between 408 MHz and 23 GHz than between 
1420 MHz and 23 GHz. For the Galactic diff'use synchrotron 
emission we would expect instead yS(o.408-23)GHz > Ai.4-23)GHz 
because of the possible steepening of the cosmic ray energy 
spectrum. In order to obtain a more quantitative estimate of the 
mean spectral index between the lower frequencies and 23 GHz, 
we used the mean value of the APS of the source- subtracted 
maps for £ e [20, 40] . In fact, in this range the APS is domi- 
nated by the diff'use synchrotron emission, whereas at higher 
multipoles the APS could still be influenced by the contribution 
of unsubtracted sources. Figure [15] shows the mean APS as a 
function of frequency for the cuts with bcut = ±40°. The spectral 
index j3 is obtained as 

< CfiVi) >^6[20,40]=< Q(V2) >^G[20,40] (Vl/V2)^^ • 
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Fig. 14. Angular power spectra of the DS-subtracted radio maps 
(green 408 MHz, fuxia 1420 MHz) extrapolated to 23 
GHz for a direct comparison with those of the WMAP 3-yr 
DS-subtracted synchrotron component (blue). The radio angu- 
lar power spectra have been smoothed to 1° to match the an- 
gular resolution of the 23 GHz map. The frequency spectral in- 
dices adopted in the extrapolation are yS(o.408-23)GHz = -2.95 and 
y^(i.4-23)GHz = -2.90 for the northern cut and y5(o.408-23)GHz = 
-2.90 and/5(i.4-23)GHz = -2.83 for the southern one. 



For the northern cut we found yS(o.408-23)GHz ~ -2.92 and 

Al.42-23)GHz ~ -2.85. For the southern cuts, yS(0.408-23)GHz ~ 

-2.76 and yS(i.42-23)GHz 2.59. These values of JS have a typ- 
ical uncertainty of a few percent according to the choice of the 
multipole range adopted to compute the mean value of the APS. 
However, it turns out that yS(o.408-23)GHz < Ai.42-23)GHz for all 
reasonable choices of the multipole range. One possible expla- 
nation for such a result is that the 23 GHz map includes one 
or more astrophysical components beside synchrotron emission. 
On the one hand, the 23 GHz map could still include some 
Galactic free-free emission, residual from the component sep- 
aration. Another and likely more relevant candidate is Galactic 
spinning dust emission. A rough estimate of the excess signal 
in the 23 GHz map can be obtained by extrapolating the radio 
results to 23 GHz, as shown in Fig. \T5\ The observed and ex- 
trapolated mean APS differ by a factor of ~ 2.5 in the northern 
cut and ~ 8.4 in the southern cut. Thus, they differ by factors 
of ~ 1.6 and ~ 2.9 respectively in terms of signal in the map. 
We repeated this calculation by computing the mean APS over 
other reasonable intervals of multipoles and found in this way 
that the uncertainty for the given values of the excess signal is 
about 20%. The difference between the observed and extrapo- 
lated values of the mean APS is smaller in the northern hemi- 
sphere, which is likely due to the compensatory contribution of 
the NPS. 

Finally, we find that the angular power spectra of the north- 
ern and southern sky are almost superimposed at 23 GHz for 
\i^cut\ ^ 40°. This can be interpreted as the combination of two 
effects. The synchrotron emission of t he NPS has a steeper fre- 
quency spectrum than the average one (lReich&Reichll988h and 
the contribution of emission processes other than synchrotron 
may be significant at 23 GHz. Therefore, the relative importance 
of the NPS with respect to the overall diffuse background dimin- 
ishes from 408 MHz and 1420 MHz to microwave frequencies. 
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Fig. 15. Mean angular power spectra, < Q(y) >^g[20,40] , of high 
latitude cuts (bcut = ±40°) against frequency. The best fit of the 
mean APS of the two lower frequencies is also plotted. 
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Table 3. Frequency spectral index of the Galactic synchrotron 
fluctuations between 408 MHz and 1420 MHz derived from the 
APS. 



4.2. Synchrotron contamination of the CMS anisotropies 

It is a standard practice to estimate the foreground contamina- 
tion of CMB anisotropies by means of the corresponding APS, 
which is usually extrapolated from the frequency range where 
the foreground component is best observed. For the Galactic 
synchrotron emission a constant spectral index in the interval 
~ [-2.5, -3.0] is commonly adopted, as suggested by the spec- 
tral behaviour of the Galactic diffuse emission at radio frequen- 
cies. We instead derive the spectral index directly from the re- 
sults of our APS analysis, thus identifying a proper value for 
each cut considered. As in the previous section, we compute the 
mean value of the APS at 408 MHz and at 1420 MHz over the 
lower multipoles (£ e [20,40]) and perform a linear extrapola- 
tion based on the two points. We prefer to work with the APS 
of the source- subtracted map rather than with the synchrotron 
power law derived by fitting it, since the former provides us with 
a value that exclusively depends on observed data. Our results 
are reported in Table [3 The obtained spectral indices vary by 
a few percent for a different choice of the multipole range used 
to calculate the mean APS. Figure [161 shows the CMB APS re- 
covered by WMAP ( Hinshaw et al. 2007), together with the syn- 
chrotron APS derived from the 1420 MHz survey, extrapolated 
to 30 GHz and 70 GHz (corresponding to the lowest and the 
highest Planck-LFI channels). We display the results obtained 
for four coverage cases (bcut = ±5°, ±20°). For comparison, we 
also extrapolated as above the APS directly extracted from the 
map for the region at |Z7| > 5° and for the all- sky. The foreground 
dominates over the CMB at 30 GHz for a wide multipole range 
if a mask excluding the Galactic plane is not applied. The fre- 
quency spectrum of free-free emission, relevant at low latitudes, 
is flatter than that of the synchrotron emission. Thus, the extrap- 
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Fig. 16. Comparison between the CMB APS retrieved bylHinshaw et alJ (|2007|) and the synchrotron angular power spectra derived 
for some cuts of the 1420 MHz map (from Table[T]and Fig.[TT]), extrapolated to 30 GHz and 70 GHz with a spectral index of -2.9 . 
Color legend (see online version): blue (upper straight lines) (northern cut), red (lower) (southern cut). The left panels 

also display the APS for the all- sky (black line) extrapolated as above. T he empty square in the top right panel marks the upper limit 
on synchrotron contamination inferred from COBE-DMR observations (iKogut et al.lll996l) . 



olated APS provides a lower limit to the overall Galactic fore- 
ground, even neglecting dust emission. 

The analysis of the low-frequency maps shows that the APS 
amplitude of the northern Galactic hemisphere is strongly in- 
fluenced by the presence of the NPS. Consequently, the results 
obtained at 1420 MHz for the northern cuts constitute a con- 
servative upper limit for the Galactic difl'use synchrotron emis- 
sion and can be used together with those of the southern cuts to 
bracket the synchrotron APS at microwave frequencies. 

At 30 GHz, a severe contamination is expected from the syn- 
chrotron emission up to ^ ~ 50 for an almost complete sky cov- 
erage (bcut = 5°). A mask excluding the region with \bgai\ < 20° 
reduces the expected synchrotron signal to about half of the 
CMB anisotropics fo r ^ > 10, whereas for lower multipoles the 
two are comparable. i Kogut et aP (11996 ) examined the COBE- 
DMR results at 31.5 GHz for \bgai\ > 20° and derived an up- 
per limit of ~ 1 lyuK on the temperature fluctuations due to syn- 
chrotron emission on angular scales of ~ 7°. This value, marked 
in Figure [161 by an empty square, is in good agreement with the 
extrapolated APS for the northern cut at 20°. 

At 70 GHz, which is the most promising channel for CMB 
anisotropy measurements since the overall foreground emission 
reaches a minimum for y ~ [60, 80] GHz (Bennett et al. 2003), 
the contribution of the Galactic synchrotron emission to the mi- 
crowave sky fluctuation field is small over the multipole range 
explored in our analysis (£ > 10). The CMB anisotropics are 
larger than the foreseen foreground fluctuations by a factor > 10 
for a cut at 5°. For bcut ~ 20° the foreground signal further 
decreases by a factor ~ 2. The extrapolation of our results to 
^ < 10 indicates that the cosmological signal should be a fac- 
tor > 2 larger than the foreground at the largest angular scales. 
The precise recovery of the CMB APS for ^ < 10 therefore re- 



mains a delicate issue, since the foreground emission is a com- 
petitive signal. However, we note that the APS extracted directly 
from the map shows a certain flattening toward lowest multi- 
poles, slightly improving the situation with respect to the above 
power law extrapolation. 

5. The APS dependence on sky position: the local 
analysis 

We have also carried out the analysis of the APS on patches of 
roughly 14?7 x 14?7, in order to describe the local variations of 
the Galactic emission at 408 MHz and 1420 MHz. Significant 
changes in the amplitude of the synchrotron APS with the con- 
sidered portion of the sky are expected, since the diff'use ra- 
dio background gradually increases toward the Galactic plane, 
where it reaches maximum intensity. 

These patches correspond to the pixels of an HEALPix map at 
^side = 4 and allow the study of the angular power spectra on the 
multipoles range ~ [60, 200-300]. An angular size Opatch ~ 14°7 
is a good compromise between the wish to divide the sky in a 
large number of areas and the need to preserve a relatively wide 
interval of statistically relevant multipoles ~ 180°/^). We 
have computed the patch angular power spectra for all the ver- 
sions of the radio maps (original, DS-subtracted and DSs only), 
both by using the HEALPix facility Anafast and by integrat- 
ing the two point correlation function (see Appendix D of La 
Porta 2007 for details). Despite the diff'erences found in individ- 
ual cases, on the average there is a good agreement between the 
angular power spectra derived with the two methods (in Fig.fTTl 
some examples of bad, fair and good cases are shown for the 
map at 1420 MHz after DS subtraction). 

The angular power spectra obtained by using Anafast typically 
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Fig. 17. Comparison between the angular power spectra de- 
rived using Anaf ast (fuxia) and via integration of the correla- 
tion function. Some examples of good (left panel), fair and poor 
(right) agreement are shown for the map at 1420 MHz after DSs 
subtraction. 



Table 4. Characteristics of the synchrotron APS best fit parame- 
ters derived in the local analysis of the 408 MHz map. refers 
to the patches covering about the brightest half of the sky, which 
includes the Galactic plane and the NPS. k^^Q correspond to the 
other half with weak high latitude emission. 



(§) 4 




parameter, x 


^min 


^max 


< x> 


O-x 


% {X^< X> +Cr^) 


a 


-3.50 


-0.70 


-2.70 


0.60 


54 


log(fe]«,/mK2) 


2.90 


5.90 


3.80 


0.81 


70 




1.00 


2.90 


2.40 


0.36 


69 



og K^oo (mK^) @ 1 420 MHz 

Fig. 20. Correlation between the best fit values obtained for the 
synchrotron emission APS normalized amplitude (^loo) at 408 
MHz and 1420 MHz. 



parameter, x 


^min 


^max 


<x> 


o-x 


% {x^< x> ±cr^) 


log (cfVmK") 


-2.30 


-0.75 


-1.50 


0.42 


66 


log (^'"VmK') 


-2.00 


0.18 


-0.86 


0.53 


72 



Table 6. Characteristics of the source term derived in the local 



analysis at 1420 MHz for \bg 

(cry 



> 45" (cD and \b,J < 45° 



parameter, x Xmin Xmax < x > cTx %(x g< x > ±o-x) 

a ^4:00 Too T80 060 70 

log(y^Ioo/mK2) -0.30 3.00 0.66 0.91 78 

log(^ino/mK^) -1-90 -0.31 -0.85 0.38 67 



The most striking result is that at each frequency the maps 
of the corresponding parameters, derived by adopting the three 
diff'erent fitting models, show a very similar morphology. Such a 
resemblance proves that the parameter patterns revealed by the 
local analysis are reliable, despite the uncertainties in the ob- 
tained parameter values. 

The slope of the synchrotron APS does not show a system- 
atic dependence on Galactic latitude, in agreement with the 
findings of Sect. 3.3. The normalized amplitude of the syn- 
chrotron APS, ^100, peaks close to the Galactic plane, which 
reflects the observed morphology. A good correlation is found 
between the normalized amplitude of the synchrotron APS at 
408 MHz and 1420 MHz, which is defined by logikf^l/mK^) ~ 
A-hBlog(y^}420/mK2), whereA = 3.15+0.02 andB = 0.88+0.02 
(see Fig.EUJ. We note that 10^ ~ (408/1420)^^ with /3 ~ -2.9, 
in agreement with the results of Table [3] 

The contribution of sources reaches a maximum in the vicin- 
ity of the Galactic plane, mainly because a less complete source 
subtraction was possible for \bgai\ < 45° than at higher latitudes 
(see Sect. 13.11) . The obtained source terms are in fair agreement 
with the values estimated by using source counts. This compari- 
son is particularly significant at 1420 MHz, where such estimates 
are more reliable. We summarize the results obtained in this case 
in Table [6l 

6. Summary and conclusions 

The aim of our analysis is to improve our understanding of the 
Galactic synchrotron emission as a foreground for CMB dedi- 
cated experiments. For this purpose, we carried out an unprece- 
dented detailed study of the Galactic radio emission, in terms of 
its angular power spectrum (APS), using total intensity all- sky 
maps at 408 MHz and 1420 MHz. 

An accurate modeling of the synchrotron APS is missing in 



Table 5. As in Tableffl but at 1420 MHz. 



present more oscillation^] and tend to be slightly flatter than the 
correlation function angular power spectra at ^ > 200. However, 
the correlation function results are less reliable at higher multi- 
poles, where the choice of the window function might have a non 
negligible influence. Consequently, we exploited the Anaf ast 
angular power spectra in the following analysis. 

5.1. Results 

The patch angular power spectra for the map after DS subtrac- 
tion are fitted exactly as done in the case of the Galactic cuts (see 
Sect. 13.3b . The maps of the obtained parameters are shown in 
Figs. [18] and [191 The results derived by using the best model are 
summarized in Tables |4]and[5l The estimated relative error of the 
synchrotron APS slope averaged over the ensamble of patches is 
\Aa/a\ ~ 25% at 408 MHz and 22% at 1420 MHz. The mean 
relative error of the normalized amplitude, kioo, is ~ 25% at 
408 MHz and ~ 20% at 1420 MHz. 



^ Anaf ast computes the APS in the Fourier space by expressing the 
temperature fluctuation field in spherical harmonics. The APS is ob- 
tained working over the whole sky, even if the map is zero outside the 
patch taken into account (the derived APS is then renormalized to the 
case of a full sky coverage). This operation is heuristically equivalent 
to computing the Fouri er transform of a disco ntinuos function, thus im- 
plying a Gibbs effect (.Arfken & Webedl2001h . 
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log (c'" / mK^) , Tl @ 1420 MHz 
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Fit no Noise 




Alpha ^„ „„i^^ , Tl @ 1420 MHz 
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Fig. 18. Maps of the best fit parameters obtained by fitting the angular power spectra of the local analysis patches at 1420 MHz. The 
maps aligned along each row refer to the same parameter. From the top, logC^^Q^^^/mK^), a^y^^^ , log(c^^^/mK^) and log(c"^'^^/mK^). 
The first and third columns correspond to the extreme cases, assuming respectively that the source contribution or the noise con- 
tamination is negligible. 



the literature so far, but is urgently required for a more precise 
and complete exploitation of the information awaited from the 
Planck satellite. It constitutes a precious input for component 
separation activities, both for the realization of spatial templates 
of the foreground and for the definition of priors on its spatial 
and frequency dependence. 

1. Being interested in the diff'use component of the synchrotron 
emission, the brighter discrete sources (DS) have been elimi- 
nated from the radio maps by 2-dimensional Gaussian fitting. 
This approach is very flexible and also permits the removal 
of extended structures. 

2. The APS was computed for both large areas and small 
patches and several consistency tests were used to check the 



reliability of the recovered APS in the case of limited sky 
coverage (see Appendix B of La Porta 2007). 
The study of the APS for various cuts, i.e. of regions with 
Galactic latitudes above or below a certain value hcut (I^cm?I ^ 
[5°, 60°]), allowed us to explore a possible dependence of the 
mean properties of the Galactic synchrotron emission on lat- 
itude, preserving at the same time the largest possible cov- 
erage, which is important when estimating the CMB APS 
because of the sampling variance. Such cuts provided infor- 
mation for t G ^max], whcrc 10 < £min ^ 30 for increas- 
ing bcut and ^^ax ~ 200 - 300 at 408 MHz and 1420 MHz, 
respectively. 

The patches correspond to the pixels of an HEALPix map at 
^side = 4, which have an angular dimension of ~ 15° and per- 
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Alpha„„^^^ , TI @ 408 MHz 



Alpha, TI ® 408 MHz 



log (c / mK^), TI @ 408 MHz 



I (c "° / mK") , TI @ 408 MHz 



log {c "° / mK"). TI @ 408 MHz 




Fig. 19. As in Fig. [HI but at 408 MHz. 



mit us to investigate the local variations of the synchrotron 
APS for multipoles larger than ~ 60. 

The derived angular power spectra were modelled in both 
cases according to Eq.[T]and a specific method was set up to 
find the best least square fit on adaptive grids of the param- 
eter space and to evaluate the uncertanties on the retrieved 
parameters (see Appendix C of La Porta 2007 for details). 
An indirect cross-check of the fit result reliability was pro- 
vided in the case of the Galactic cuts by the estimated source 
terms, which are consistent with the expectations from extra- 
galactic source counts at both frequencies. Nowadays, source 
counts at 1 .4 GHz are well established down to very low flux 
limits. It is remarkable that the source angular power spectra 
obtained independently from source counts and from the fit 
of the survey APS are in good agreement. 
The slope of the synchrotron APS, Q ~ H^, changes with 
bcut without showing a well-defined regular trend, although 



Fit no Noise 



mK^), TI @ 408 MHz 




Alpha , TI @ 408 MHz 



'm 



log (c "° / mK^) , TI @ 408 MHz 




it is found to be typically steeper for hgai ^ 20°. For the 
cuts, a varies in the range ~ [-3.0, -2.6] at both frequen- 
cies. However, the analysis of the small patches gives ev- 
idence that locally the synchrotron APS can be much flat- 
ter for t > 60, reaching in some cases values of a 0.8. 

The normalized amplitude, ^loo = k x 100^, gradually in- 
creases toward the Galactic plane, following, as expected, 
the background radio emission gradient. A good correlation 
exists between the results obtained for ^loo at 408 MHz and 
1420 MHz, for both cuts and patches. This is expected, given 
that the spectral properties of the electron density distribu- 
tion responsible for the Galactic difl'use non-thermal emis- 
sion should be the same in that frequency range and further 
supports the reliability of the obtained estimates of the syn- 
chrotron APS. 

5. The maps of ^loo and a resulting from the local analysis rep- 
resent the starting point for the simulation of small-scale 
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fluctuation fields to be added to the DS-subtracted maps 
to build phenomenological templates of the Galactic syn- 
chrotron emission. At present, an empirical approach is the 
most reliable way to proceed in the realization of realistic 
templates of the foreground, given the poor knowledge of the 
Galactic magnetic field and of the cosmic ray electron den- 
sity distribution needed for a 3 -dimensional physical mod- 
elling. 

The issues (discussed in Sect. 14.11) related to the K-ban d syn- 
chrotron component retrieved by iHinshaw et al.l (|2007|) are a 
clear example of the difficulties encountered in the foreground 
separation due to the lack of adequate priors and guess. We per- 
formed a source subtraction on that map and produced a map 
of the Galactic diffuse non-thermal emission at 23 GHz to be 
compared with those at lower frequencies. 

6. The extrapolation to 23 GHz of the APS obtained at 408 
MHz and 1420 MHz for higher latitude regions (\bgai\ > 
40°) reveals that the mean spectral index (Q(y) oc y^^) 
Ao.408-23)GHz ^ Ai.4-23)GHz, which is the oppositc of what is 
expected for synchrotron emission. We estimate the excess 
of the signal in the 23 GHz map to be > 50% by using the 
mean value of the APS at lower multipoles. This result can 
be interpreted in terms of additional contributions to the 23 
GHz Galactic non-thermal emission, which could be mainly 
due to anomalous dust. 

A direct application of the presented analysis is to determine 
the level of contamination of the CMB anisotropics due to the 
Galactic diffuse synchrotron emission at different angular scales. 
The conclusions reported below refer to the cuts, which are more 
relevant for CMB measurements because of the large coverage. 
However, similar considerations could be repeated using the re- 
sults of the local analysis, thus allowing the identification of the 
clearest sky areas, which is essential for the success of ground- 
based CMB experiments. 

7. An important - although not unexpected - outcome of the cut 
analysis is that the amplitude of the APS for the northern 
hemisphere cuts above ~ 20° is raised by the presence of the 
NPS. Consequently, the results obtained separately for the 
two Galactic hemispheres were used to bracket the APS of 
the synchrotron emission. 

8. We used the APS results at 408 MHz and 1420 MHz to de- 
termine the frequency spectral index to be adopted in the ex- 
trapolation to the microwave range for each coverage case. 
In particular, we found that /5(o.408-i.4)GHz ^ [-3.2, -2.9] with 
an uncertainty of a few percent. 

9. The extrapolation to 30 GHz of the synchrotron APS 
obtained at 1420 MHz led to a signal in good agreement 
with the upper limit fixed by COBE-DMR, thus supporting 
the reliability of our results. At this frequency the syn- 
chrotron emission constitutes a severe contamination of 
CMB anisotropics at the largest angular scales (£ < 40). 
Nevertheless, a cut at ~ 20° would reduce the synchrotron 
fluctuations to about half the cosmological ones. The same 
holds at y ~ 70 GHz for a cut at ~ 5° and the situation 
further improves when excluding a larger portion of the sky 
around the Galactic plane. This implies that even though 
the current treatment of the foregrounds does not permit an 
accurate removal of the synchrotron emission, the latter does 
not prevent the recovery of the bulk of the cosmological 
information encoded in the CMB temperature APS. 



A deeper understanding of the foreground is indispensable to 
settle other important issues in the perspective of the forthcom- 
ing Planck mission. Planck should achieve a sensitivity com- 
parable to the cosmic variance and its performance should be 
limited mainly by foregrounds, thanks to the e xtremely accu- 
rate control of all in strumental systematic effects (iBurigana et al.l 
120041: iMennefla et al.ll2004^) . 

A first issue is the estimate of the CMB temperature-polarization 
APS, which carries information about the reionization history 
and the tensor- to- scalar ratio (see, e.g. Kogut 2003, Kogut et 
al. 2003, and references therein). Possible foreground residual 
contamination in the total intensity CMB anisotropy map would 
affect fine analysis based on the estimate of the cross-correlation 
APS, also because the polarized component of the cosmological 
signal is orders of magnitude lower. 

Another issue is the evaluation of the Gaussianity of the 
primordial fluctuations, which in the standard inflationary 
paradignfl generate the structures observed in the Universe to- 
day. Gaussianity tests are a powerful tool, complementary to 
the tests exploiting the APS, which allow us to probe the "con- 
cordance" model ( Spergel et al. 20 07b and also to distinguish 
among inflationary models (IBartolo et al.l l2004l) . The level of 
non-Gaussianity predicted by the "concordance" model can- 
not be detected by WMAP, but in principle it should be ob- 
servable with Planck. Galactic foregrounds are non-Gaussian 
and anisotropic and even low-level contamination in the maps 
can produce detectable non-Gaussianities (see, e.g. Naselsky et 
al. 20 05), although they ha ve minimal effects on the estimated 
APS (iHinshaw et al.l l2003l) . Consequently, the foreground re- 
moval has to be extremely accurate, so as not to limit Planck 
in verifying this cruciaQ prediction. 

Acknowledgements. We are grateful to R. Wielebinski for a careful reading of 
the original manuscript. We wish to thank G. De Zotti, L. Toffolatti and R. 
Rebolo for helpful discussions. L.L.P. warmly thanks R. Wielebinski and A. 
Zensus for granting a post-doc fellowship. C.B. acknowledge the support by the 
AST contract "Planck LFI Activity of Phase E2". We are grateful to M. Genghini 
for technical supp ort. Some of the re sults in this paper have been derived us- 
ing the HEALPix dGorski et al.1l2005l) package. The availability of the WMAP 
3-yr maps is acknowledged. We warmly thank the anonymous referee for useful 
comments. 



References 

Alishahiha, M., Silverstein, E., & Tong, D. 2004, Phys. Rev., 70, 12 

Arfken, G. B., & Weber, H. J., 2001, Mathematical Methods for Physicists, 

Harcourt Academy Press, p. 893 
Arkani-Hamed, N., Creminelli, P., Mukohyama, S., & Zaldarriaga, M. 2004, 

JCAP, 4, 1 

Baccigalupi, C, Burigana, C, Perrotta, R, et al. 2001, A&A, 372, 8 

Banday, A. J., & Wolfendale, A. W. 1990, MNRAS, 245, 182 

Banday, A. J., & Wolfendale, A. W. 1991, MNRAS, 252, 462 

Bartolo, N., Komatsu, E., Matarrese, S., & Riotto, A. 2004, Phys. Rev., 402, 103 

Bennett, C. L., Hill, R. S., Hinshaw, G., et al. 2003, ApJS, 148, 97 

Bouchet, F. R., Prunet, S., & Sethi, S. 1999, MNRAS, 302, 663 

Bouchet, F. R., & Gispert, R. 1999, NewA, 4, 443 

Burigana, C, Finelli, F., Salvaterra, R., et al. 2004, MmSAI, 5,415 

Burigana, C, La Porta, L., Reich, R, & Reich, W. 2006, AN, 327, 491 



^ Simple (standard) inflationary scenarios (see the reviews by Lyth 
& Riotto 1999 and Linde 2005) predict the existence of Gaussian den- 
sity fluctuations and of gravitational waves with a nearly scale-invariant 
spectrum. 

^ Standard cosmologies predict a minimum level of non-Gaussianity 
f or the the pri n iordial perturbations. In alternative cosmologies 
dLvth et all l2003l : lArkani-Hamed et all l2004l : lAlishahiha eFaP l2004h 
such a lower limit is even higher. Consequently, detection or non- 
detection of non-Gaussianities sheds light on the physics of the early 
Universe. 



La Porta et aL: APS analysis of total intensity all-sky surveys 



15 



Chepurnov, A. V. 1998, A&AT, 17, 281 

Chiang, L.-Y., Coles, P., & Naselsky, R D. 2007, JCAP, 1, 21 

Cho, J., & Lazarian, A. 2002, ApJL, 575, 63 

Cho, J., & Lazarian, A. 2003, NewAR, 47, 1143 

Davies, R. D., Dickinson, C, Banday, A. J., et al. 2006, MNRAS, 370, 1125 
De Oliveira-Costa, A, Tegmark, M., Davies, R. D., et al. 2004, ApJL, 606, 89 
Dickinson, C, Davies, R. D., & Davis, R. J. 2003, MNRAS, 341, 369 
Draine, B. T., & Lazarian, A. 1998, ApJL, 494, 19 
Egger, R. J., & Aschenbach, B. 1995, A&A, 294, L25 
Franceschini, A., Toffolatti, L., Danese, L., De Zotti, G. 1989, ApJ, 344, 35 
Finkbeiner, D. R, Davis, M., & Schlegel, D. J. 1999, ApJ, 524, 867 
Finkbeiner, D. R 2003, ApJS, 146, 407 
Finkbeiner, D. 2004, ApJ, 614, 186 

Giardino, G., Banday, A. J., Gorski, K. M., et al. 2001, A&A, 371, 708 
Gorski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 
Haslam, C. G. T. 1974, A&AS, 15,333 

Haslam, C. G. T., Salter, C. J., StofFel, H., & Wilson, W. E. 1982, A&AS, 47, 1 
Hildebrandt, S. R., Rebolo, R., Rubino-Martin, J. A., et al. 2007, MNRAS, sub- 
mitted, astro-ph/07061873 
Hinshaw, G., Spergel, D.N., Verde, L. et al. 2003, ApJS, 148, 135 
Hinshaw, G., Nolta, M. R., Bennett, C. L., et al. 2007, ApJS, 170, 288 
James, F., & Roos, M. 1975, Comput. Phys. Comm., 10, 343 
Jamrozy, M. 2004, A&A, 419, 63 

Kogut, A., Banday, A. J., Bennett, C. L., et al. 1996, ApJL, 464, 5 
Kogut, A. 2003, NewAR, 47, 977 

Kogut, A., Spergel, D.N., Barnes, C. et al. 2003, ApJS, 148, 161 

La Porta, L., Reich, P, Burigana, C, & Reich, W. 2005, MPIfR-Memo,l 

La Porta, L. 2007, PhD thesis, Bonn University, Germany 

Linde, A. 2005, NewAR, 49, 35 

Lyth, D. H., & Riotto, A. A. 1999, Phys. Rep., 314, 1 

Lyth, D. H., Ungarelli, C, & Wands, D. 2003, Phys. Rev., 67, 2 

Mennella, A., Baccigalupi, C, Balbi, A., et al. 2004, Recent Research 

Developments in Astronomy & Astrophysics, Vol. 2, astro-ph/0402528 
Naselsky, P D., Chiang, L.-Y., Olesen, P, & Novikov, I. D. 2005, Phys. Rev., 

72, 063512 

Naselsky, P D., Novikov, I. D., & Chiang, L.-Y. 2006, ApJ, 642, 617 
Naselsky, P D., & Verkhodanov, O. V. 2006, MNRAS, submitted, 
[astro-ph/0609409 

Negrello, M., Gonzalez-Nuevo, J., Magliocchetti, M., et al. 2005, MNRAS, 358, 
896 

Negrello, M., Magliocchetti, M., & De Zotti, G. 2006, MNRAS, 368, 935 

Paladini, R., De Zotti, G., Davies, R. D., & Giard, M. 2005, MNRAS, 360, 1545 

Platania, P, Burigana, C, Maino, D., et al. 2003, A&A, 410, 847 

Prandoni, I., Gregorini, L., Parma, P. et al. 2001, A&A, 365, 392 

Reich, P, & Reich, W. 1986, A&AS, 63, 205 

Reich, P, & Reich, W. 1988, A&AS, 74, 7 

Reich, P, Testori, J. C, & Reich, W. 2001, A&A, 376, 861 

Reich, P, Reich, W., & Testori, J. C. 2004, Proceeding The Magnetized 

Interstellar Medium, Antalya, Turkey, Eds: B. Uyaniker, W. Reich, and R. 

Wielebinski, Copernicus GmbH, Katlenburg-Lindau, p. 63, 2004 
Reich, W. 1982, A&AS, 48, 219 

Rybicki, G. B., & Lightman, A. P. 1979, Radiative processes in astrophysics. 

New York, Wiley-Interscience 
Salter, C. J. 1983, BASI, 11, 1 

Spergel, D. N., Bean, R., Dore, O., et al. 2007, ApJS, 170, 377 

Strong, A. W., Moskalenko, I.V, & Ptuskin, V. S. 2007, Annu. Rev. Nucl. Part. 

Sci., 57, 285 
Tauber, J. 2004, AdSpR, 34, 91 

Tegmark, M., & Efstathiou, G. 1996, MNRAS, 281, 1297 
Tegmark, M., Eisenstein, D. J., Hu, W, & De Oliveira-Costa, A. 2000, ApJ, 530, 
133 

Toffolatti, L., Argueso Gomez, R, De Zotti, G. et al. 1998, MNRAS, 297, 117 
Testori, J. C, Reich, P, Bava, J. A. et al. 2001, A&A, 368, 1123 



